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Abstract 



A non perturbative finite size scaling technique is used to study a 
running coupling in lattice Yang-Mills theory coupled to a bosonic Wilson 
^~») , spinor field in the Schrodinger functional scheme. This corresponds to 

f"^ ■ two negative flavours. The scaling behaviour in this case is compared to 

f-^. ' quenched results and to QCD with two flavours. The continuum limit is 

f*"*) , confronted with renormalized perturbation theory. 

On 
On 



1 Introduction and motivation 



One of the main goals of the ALPHA collaboration is the non perturbative com- 
putation of the strong coupling as for energy scales ranging from hadronic scales 
fH \ to perturbative high energy scales. To this end a finite volume scale dependent 

renormalization scheme for QCD has been invented [jij g]. Its central object is 
the step scaling function £ which describes the running of the coupling under 
■ a discrete change of the scale. Other ingredients include Schrodinger functional 

boundary conditions, 0(a) improvement and non perturbative renormalization. 
This method has been used successfully in the quenched approximation (see || 
for a review). First results for full QCD with two flavours have been obtained 
recently Q . In figure [I] the results for the step scaling function at the coupling 
g = 0.9793 are shown and compared with the quenched approximation and the 
perturbative continuum limitjj]. Here an O(a) improved action was used. The 
observed cut-off effects are larger than in the quenched approximation. A naive 
extrapolation of the Monte Carlo data to the continuum limit yields a value 
which lies 2-3% above the perturbative estimate, although the coupling g 2 is 
expected to be small enough to be in the perturbative regime. Furthermore, 

"Talk given at the International Symposium on Lattice Field Theory, June 28-July 3, 1999, 
in Pisa, Italy. 

lr rhese data were obtained in collaboration with A. Bode, R. Frezzotti, M. Guagnelli, K. 
Jansen, M. Hasenbusch, J. Heitger, R. Sommer and P. Weisz. 
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Figure 1: Preliminary results for the step scaling function in QCD with two 
Eavours at u — 0.9793 compared with the quenched approximation. Through- 
out this paper we denote Nf — 2 by circles, Nf = by squares and Nf = —2 
by triangles. 

the cut-off effects in the quenched and in the full theory computed in lattice 
perturbation theory to two loops are of the same size These observations 
can be interpreted as follows: 

• statistical fluctuation 

• true 3% deviation from renormalized perturbation theory 

• lattice artifacts too large for continuum extrapolation 

The investigation of these possibilities in full QCD would consume a substantial 
amount of computer time. Therefore we use Yang-Mills theory coupled to a 
bosonic Wilson spinor field, which corresponds to setting Nf = —2 in the QCD 
partition function, as a toy model to study the extrapolation to the continuum 
limit for a system different from pure gauge theory. This model has been called 
the bermion model in the literature [|| Q . 
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2 The bermion model in the Schrodinger func- 
tional setup 

Let the space time be a hypercubic Euclidean lattice with lattice spacing a and 
volume L 3 x T. In the following we set T = L. The SU(3) gauge field U(x, fi) is 
defined on the links while the bermion field 4>{x) which is a bosonic spinor field 
with color and Dirac indices is defined on the sites of the lattice. In the space 
directions we impose periodic boundary conditions while in the time direction 
we use Dirichlet boundary conditions. The boundary gauge fields can be chosen 
such that a constant color electric background field is enforced on the system 
which can be varied by a dimensionless parameter 77 JlJ . 

As explained above our goal is to continue the exponent of the fermion 
determinant in the QCD partition function to the negative value Nf = —2. This 
can be achieved by integrating the bosonic field <fi with the Gaussian action 

Sb = ^ |M0(V)| 2 , where (I) 

X 

^-M^x) = (V + m )4,(x). (2) 

IK 

T> is the Wilson Dirac operator with hopping parameter k = (8 + 2mo) . For 
the gauge fields we employ the action 

S G = ^52w(pMl-U(p)). (3) 
9q p 

The weights w (p) are defined to be one for plaquettes p in the interior and they 
equal q for time like plaquettes attached to the boundary. The choice c t = 1 
corresponds to the standard Wilson action. However, c t can be tuned in order 
to reduce lattice artifacts. For the bermion case we have in this work always 
chosen c t = I . 

Now the Schrodinger functional is defined as the partition function in the 
above setup: 

Z = J VUV<j)V(j) + e - SG - SB (4) 
= JvU e ~ Sa ^ det{M+M) N f/ 2 , (5) 
with Nf = —2. The effective action 

r = -io g z (6) 

with the perturbative expansion 

r = . 9 - 2 r + r a + , 9o 2 r 2 + . . . (7) 
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is renormalizable with no extra counterterms up to an additive divergent con- 
stant. That means that the derivative T' = ^ is a renormalized quantity and 

.9 2 (L) = ryrv (8) 

defines a renormalized coupling which depends only on L and the bermion mass 
m. Note that this coupling can be computed efficiently as the expectation value 

7fn = (if)' ^ e mass m ^ S defined via the PC AC relation ||. Here we use 
the fermionic boundary states of the Schrodinger functional^ to transform this 
operator relation to an identity (up to 0(a)) for fermionic correlation functions 
which can be computed on the lattice. This gives a time dependent mass m(xo)- 
The mass m is then defined by 

m(f ) T even, , , 

i ( m (2^i) +m (Z+i)) Todd. K ) 

To define the step scaling function cr(s.u) let u — g 2 {L) and m(L) = 0. 
Then we change the length scale by a factor s and compute the new coupling 



as 



" — 9 2 { S L)- The lattice step scaling function £ at the resolution L/a is defined 



S(*,u,o/L) = f{sL)\_ HL)=Ui m(L)=0 . (10) 

Note that the two conditions fix the two bare parameters go and k. The con- 
tinuum limit cr(s,u) can be found by an extrapolation in a/L. That means the 
computational strategy is as follows: 

1. Choose a lattice with L/a points in each direction. 

2. Tune the bare parameters go and k such that the renormalized coupling 
g 2 {L) has the value u and m(L) = 0. 

3. At the same value of go and k simulate a lattice with twice the linear size 
and compute u' — g 2 (2L). This gives E(2,u, -g). 

4. Repeat steps 1.-3. with different resolutions L/a and extrapolate f- — > 0, 
which yields <r(2, u). 



3 Results 

We have performed Monte Carlo simulations on APEIOO/Quadrics parallel com- 
puters with SIMD architecture and single precision arithmetics. The size of the 
machines ranged from 8 to 256 nodes. The gauge fields and the bermion fields 
have been generated by hybrid ovcrrclaxation including microcanonical reflec- 
tion steps. While we measured the gauge observables after each update of the 

2 Fermionic observables are constructed independently of the number of dynamical flavours 
N f . 
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fields the fermionic correlation functions were determined only every 100th to 
1000th iteration, since their measurement involves the inversion of the Dirac op- 
erator. The statistical errors have been determined by a direct computation of 
the autocorrelation matrix. The largest run has been performed for the lattice 
size L = 24 which took about 12 days on the largest machine. 
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Figure 2: Results for the lattice step scaling function for Nf = — 2 at u = 
0.9793 confronted with renormalized perturbation theory. 

In figure || our results for the bermion lattice step scaling function 
E(2, 0.9793, a/L) are shown for resolutions 4, 5, 6, 8 and 12 and compared with 
renormalized perturbation theory. The lattice artifacts are consistent with O(a) 
effects. Included is a linear fit to the data which extrapolates to a continuum 
limit which lies above the perturbative estimate but is compatible with it. Fur- 
thermore a linear plus quadratic fit with the constraint that the continuum 
limit is given by the perturbative estimate is shown. All the data points are 
compatible with this fit. Note that L/a — 4 has been ignored in these fits. 

In figure ^ we compare these data with results from the quenched approxima- 
tion with the standard Wilson action and with a perturbatively O(a) improved 
action. Again the most naive extrapolation of the unimproved data would ex- 
trapolate to a value slightly above the perturbative estimate. However, if we use 
universality, i.e. the agreement of the continuum limit of the two data sets as 
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Figure 3: Results for the lattice step scaling functions for Nf — in the 
unimproved and the improved case at u — 0.9793 confronted with renormalizcd 
perturbation theory. 



a constraint, their joint continuum limit is fully compatible with perturbation 
theory. 

Furthermore we observe that in the bermion theory the cut off effects are 
larger than in the quenched approximation. 

Since the extrapolation to the continuum limit in the quenched approxima- 
tion is much easier in the improved case we plan to study the O(a) improved 
bermion model. This will also allow us to compare with the fermionic theory 
on equal footing. 

This work is part of the ALPHA collaboration research program. We thank 
DESY for allocating computer time to this project and the DFG under GK 271 
for financial support. 
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